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Fractional charge revealed in computer simulations of resonant tunneling in the 

fractional quantum Hall regime 
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The concept of fractional charge is central to the theory of the fractional quantum Hall effect 
(FQHE). Here I use exact diagonalization as well as configuration space renormalization (CSR) 
to study finite clusters which are large enough to contain two independent edges. I analyze the 
conditions of resonant tunneling between the two edges. The "computer experiment" reveals a 
periodic sequence of resonant tunneling events consistent with the experimentally observed fractional 
quantization of electric charge in units of e/3 and e/5. 
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Perhaps, the most intriguing feature of the FQHE 1 is 
the existence of quasiparticles whose electric charge is a 
simple fraction of the elementary charge e. 2 Quasiparti- 
cles of charge e* = e/3 and e/5 have been first observed 
experimentally in the v = ^ and v = | fractional states, 
respectively, using resonant tunneling via a quantum an- 
tidot (a potential hill). 3-5 

Since the bulk fractional state is an insulator, the most 
interesting transport properties of the system are asso- 
ciated with the edges, particularly with tunneling into 
or between the edges. 3-8 Due to cluster size limitations, 
computational studies of edge physics have focused on 
the properties of a single edge, such as non-universality 
of the tunneling exponent 9 ' 10 or reconstruction of the 
charge density. 11-14 The properties of the Laughlin wave- 
function describing a dual-edge system have been studied 
in cylindrical 15 and disk geometries. 16 Study of edge to 
edge tunneling through a bulk fractional state requires 
clusters large enough to contain two independent edges. 

Exact diagonalization (ED) of finite clusters has been 
very fruitful in helping to understand the physics of 
FQHE. 2 ' 9-16 Ordinary electronic structure methods fail 
for this system because the kinetic energy is quantized by 
the magnetic field. ED, or "Full CI" in quantum chemi- 
cal terminology, imposes severe restrictions on the cluster 
size, since the dimensionality of the Hilbert space grows 
exponentially with the number of particles N: 
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Here L = N/f, f < 1 is the filling factor. The quantity 
in brackets reaches maximum value of 2 at f=h- 

In the past we have perfected the Lanczos 
technique, 17,18 both Hermitian 19 ' 20 and not, 21-23 to work 
with matrices up to ~ 10 9 x 10 9 . Equation (1) translates 
this into about iV=12 particles at f=\ or 7V=16 at f=\- 
Whereas exact solutions for up to N=22 are sometimes 
possible, 20 N > 12 normally require approximate meth- 
ods. Here I use ED and also an approximate method to 
model the resonant tunneling experiments. 3 5 



In a periodic sequence of resonant tunneling events 
was observed as either the magnetic field H or the back- 
gate voltage Vbg were varied. The tunneling events are 
thought of in terms of a quasiparticle tunneling through 
the bulk of the fractional state between the outer edge 
of the sample and the inner edge formed around the an- 
tidot. The periodicities AH and AVbg were related 3 to 
the quasiparticle charge e*. 
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Fig. 1 Exact ground-state energy of N = 12 electrons in 
2D confined by Coulomb attraction to a uniformly-charged 
annulus (inset) of charge density cr = | and | of the filled 
Landau level. The units are e 2 /in- The upper panel also 
shows CSR results for K = 10 4 and 10 5 . Full Hilbert space 
dimensionality is 10 8 — 10 9 and varies with M. 26 The § data 
is shifted upward by 1.5. The lower panel exposes the FQHE- 
related structure in Em by subtracting the greatest convex 
minorant Em [the dotted purple line]. 

In order to mimic the experimental setup I consider 
a planar FQHE sample with two unconnected edges (in- 
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set in Fig. 1). A electrons in the lowest Landau level 
are confined by the potential of a uniformly-charged disk 
with a hole in the center, positioned in the plane of the 
two-dimensional (2D) electron gas. The positive charge 
density a and the inner radius R\ of the disk are free pa- 
rameters. The outer radius i?2 is always chosen such that 
the whole system is neutral. The electronic density p{r) 
confines itself between i?i and i?2, falling off sharply be- 
yond this range. Setting a to a fraction v = -|, |, etc. of 
the density o\ of the completely filled Landau level con- 
trols the fractional state, with p(r) approaching vo\ (for 
N — > oo) far from both edges. Near the edges p(r) is 
known to exhibit oscillatory behavior thought to decay 
slowly into the bulk. 13 Such behavior prevents formation 
of a well-defined density plateau between the edges in the 
finite clusters studied here numerically. 

An increase in R\ strengthens the antidot and expels 
charge from inside of the antidot towards the outer edge. 
The charge expelled does not accumulate in the bulk be- 
cause of neutrality considerations and because of incom- 
pressibility of the bulk fractional state. I prefer to use the 
"missing charge" Q — cnrR\/e as a variable, instead of 
R\. When Q is continuously increased, the ground state 
of the system reconstructs via a step-like process. The 
reconstruction events correspond to ground state degen- 
eracies, when it costs no energy to transfer charge from 
the inner to the outer edge. This is precisely the condi- 
tion for resonant tunneling through the antidot. 

In the disk geometry the single-particle states ip m in 
the lowest Landau level are characterized by the angu- 
lar momentum m = 0,1,..., and the total angular mo- 
mentum M — Y^, m is conserved. The Coulomb matrix 
elements are known. 24,25 The matrix elements of the con- 
fining potential are V m = V m (R2) — V m {Ri), where 

V m (R)= [ [ rf 2 pA^|^ m (r)| 2 (2) 
J J P <R I r -Pl 

m 

= {2^eal H Ze- z Y,[^rMZ) + q+MZ)]Z\ 

i=0 



Here £h — yhcfeH is the magnetic length, Z = 
R 2 /M% q± = 1, g±, = (2m ± l)q , m -i, qf m = (2m - 

2* ± lktm-i + 2g£i, m -i - 2gjL lim _ 1 . 

For a given set of A, Q, and a I find the lowest en- 
ergy Em(Q) at each M. The ground state energy is then 
E(Q) = mini^M(Q). The ground state reconstruction 
events occur via level crossings of branches with different 
M and lead to a step-wise function M{Q). The number 
p of the steps that occur per AQ = 1 may be related to 
the charge e/p that is moved from the inner to the outer 
edge per one reconstruction event. 

Figure 1 shows E M (Q) for A = 12, Q = 2. The se- 
quence of sharp cusps on the right curve are the ^ frac- 
tional states. The state at M — 270 is the true ground 
state at Q ~ 2 and a — |o"i, whereas the states with 
M = 270 ± 12, M = 270 ± 24, etc. are candidates for the 
ground state at different Q. The quasi-periodicity with 



AM = 12 is due to the approximate invariance of the an- 
tidot Hamiltonian with respect to the Laughlin's quasi- 
hole creation operator 2 Ao, also known to be the gener- 
ator of infinitesimal magnetic translations. 27 Applied to 
an arbitrary many-electron wave function 'F, it translates 
it in the angular momentum space, m — ► m+1. The total 
angular momentum then transforms as M — > M + N: 

*M+W « Aq^m (3) 

At Q = the ground state occurs, approximately, at the 
Laughlin's angular momentum 28 

M*(Q = 0) = N{N 2 - 1) . (4) 

I generalize this for a disk with the missing charge Q as 
= (N + Q)(N + Q-1) _ Q(Q - 1) 

For N = 12, Q = 2 I get M* /3 = 270 and M* /5 = 225 
(cf. Fig. 1). 

The single-particle orbitals ip m (r) are localized near 
r = £.H\/2m, where £h is the magnetic length. There- 
fore, in a macroscopic system whose density approaches 
a constant vai in the bulk, the operator Ao pushes the 
density out of the center, creating an effective positive 
charge e* = ve. This is an exact formal property of ^4 
but relates to the physical system via the approximate 
invariance (3). Indeed, Eq. (5) can be obtained from (4) 
by applying A^ /u) : 

M* = M*(Q = 0) + NQ/v. (6) 

Figure 1, therefore, suggests that the ground state of 
the | system changes in steps of AM = N, transfer- 
ring charge e* = e/3 from the inner to the outer edge at 
every step. According to (6) the steps should occur at 
AQ «i/ = §. 

Remarkably, the range of M that corresponds to the | 
fractional state exhibits double periodicity, AM = 6 = 
A/2. This should cause branch crossings twice as of- 
ten, AQ w vj2 = 1/5. Consequently, the charge trans- 
ferred per one reconstruction event is e* = ve/2 = e/5, 
in agreement with the experiment. 4 The double period- 
icity in Em{Q) at a — |cri shows up also for N = 11, 10, 
9, and 8, though it is less pronounced for smaller A. 

Figure 1(b) exposes the tiny structure in Em(Q) by 
subtracting its greatest convex minorant Em [dotted pur- 
ple line in Fig. 1(a)]. I notice that this structure, which is 
the manifestation of the FQHE, is insensitive to the con- 
fining potential, and is practically the same for a = |cri 
and gal, suggesting rigidity of the wave function with re- 
spect to the confining potential [8^ M/5V(r) being small 
by some measure], whereas the principal effect of V(r) 
is to select which M is the ground state. Figure 1(b) 
also suggests that traces of AM — 4 periodicity may be 
present in the area near M^ /7 = 210. By the same logic, 
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AM = N/3 at | filling corresponds to e* = e/7, although 
larger clusters are needed to separate the | and | states. 
Indeed, from Eq. (4), the condition My 5 — M*p > N 
leads, at Q = 0, to N > 13. 

I further use the CSR approach to study larger clus- 
ters up to N—15 and to compute Em(Q) for several 
values of Q and obtain the branch crossings directly. 
The CSR reduces the dimensionality of the many-body 
Hilbert space by selecting its most relevant subspace. It is 
similar in spirit to the basis set reduction ( "BSR" ) tech- 
nique by Wenzel and Wilson 29 and a number of related 
algorithms. 30-33 The CSR iteratively improves a basis set 
of K many-body configurations using their weight in the 
solution as the criterion of relevance. 

The CSR algorithm works as follows. I start with an 
arbitrary set of K many-body configurations and diago- 
nalize the Hamiltonian in the subspace they span. The 
result is a state vector expanded in many-body config- 
urations. I then retain K' < K configurations by dis- 
carding the ones which have little weight. I re-expand 
the set with the new configurations that have large ma- 
trix elements with those retained. When repeated, the 
procedure converges after 10 — 15 iterations (Fig. 2) and 
yields some optimal set of configurations. 




10 20 30 



Iteration 

Fig. 2 Typical CSR energy convergence, N = 12, a = §cri, 
M = 216, K = 10 5 . The inset shows convergence with K 
when K is increased by a factor 1.2 at every iteration. "+" 
(red) and "x" (green) data sets indicate insignificance of the 
starting value of K except for the initial iterations. The stars 
(blue) give overlap with the exact eigenvector. The squares 
(purple) show the magnitude of projection of the exact eigen- 
vector onto the CSR subspace (that is, maximum overlap with 
any vector in the subspace) , and as such characterize the qual- 
ity of the subspace. 34 

The resulting basis truncation is essentially many- 



body, and cannot be achieved by truncating or rotating 
the single-particle basis. ED performed in the subspace 
yields a variationally-stable ground state energy that con- 
verges to the exact value as K is increased. In practice, 
I do not keep K constant, but increase it from iteration 
to iteration (Fig. 2, inset), monitor the convergence, and 
extrapolate as 1/K — > 0. 35 

Figure 1 compares CSR results for K = 10 4 and 10 5 
against the exact solution. We see that the essential 
FQHE structure survives the basis truncations of sev- 
eral orders of magnitude, and that a reasonable accuracy 
is achieved for K = 10 5 in the whole range of M. Qual- 
itative results are obtained already with K as small as 
10 4 . 

I used CSR to compute Em{Q) for Q in steps of 0.25, 
interpolated between these points with a cubic spline and 
found mm Em (Q) over M at every Q (Fig. 3). I used 
K = 200,000 for N < 13, K = 500,000 for N = 14, 
and K = 1, 000, 000 for N = 15. K was doubled in some 
calculations where the extrapolation to 1 /K — » seemed 
ambiguous. 
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Fig. 3 Angular momentum M(Q) of the ground state 
changes in steps as the "missing charge" Q is tuned continu- 
ously. The steps occur at the ground state degeneracies, when 
it costs no energy to move a quasiparticle from the inner to 
the outer edge, and can be associated with the resonant tun- 
neling between the edges through the bulk of the FQH state. 

Figure 3 shows steps M(Q) for N from 11 through 
15, as the "missing charge" Q is continuously increased. 
The left panel shows that three steps typically occur per 
AQ = 1 for all N. The general slope is consistent with 
Eq. (6), and most of the steps in the left panel have 
AM = N precisely. The data for the | fraction (the 
middle panel) shows steps occurring about twice as of- 
ten. These observations are in line with our expectations 
based on the quasi-periodicities seen in Fig. 1. Care- 
ful examination of the | data shows that the change 
M — ► M + N occurs usually in two steps, although these 
steps are not always equal to N/2. The data for | has 
been computed only for N > 13 as discussed above. It 
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shows behavior similar to the | data. An expectation of 
7 steps per AQ = 1 (or AM — N/3 as we would expect 
for e* = e/7) cannot be confirmed. This could mean that 
N = 15 cluster size is not large enough to distinguish the 
| fraction. It could also indicate a genuine property of 
the | fraction that needs to be understood. 

In conclusion, I have conducted a "computer exper- 
iment" designed to model the key elements of the real 
experiment on resonant tunneling through a quantum 
antidot. 3 The data on small clusters obtained using ED 
and CSR reveal a sequence of the ground state recon- 
struction events consistent with the periodicity of the 
resonant tunneling peaks observed experimentally. 

The CSR approach employed here offers a number of 
advantages in modeling many-body clusters. In particu- 
lar, truncation of the single-particle basis, whenever pos- 
sible (such as restricting the range of m 26 ) occurs auto- 
matically within CSR, which discards irrelevant many- 
particle configurations and thus effectively removes any 
single-particle state that contributes to none of the rel- 
evant configurations. This feature itself can be a major 
simplification, because the relevance of a single-particle 
orbital cannot be always judged in advance. 

The term "renormalization" I use reflects the spirit 
of the renormalization theory, yet I have not observed 
any defined fixed point: the set of K relevant configu- 
rations keeps changing slightly upon convergence, with 
some marginally-relevant configurations being replaced 
with other similarly relevant. 

In general, performance of CSR may benefit if it is 
augmented with rotations of the single-particle basis, 36 
though the choice of the basis is fixed by symmetry in 
the case of the lowest Landau level in disk geometry. 

I would like to thank V.J. Goldman for numerous en- 
lightening conversations. I also acknowledge stimulating 
discussions with A.L. Aleiner. 
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